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Long-range properties of the two-point correlation function of the electromagnetic field produced 
by an elementary particle are investigated. Using the Schwinger-Keldysh formalism it is shown 
that this function is finite in the coincidence limit outside the region of particle localization. In this 
limit, the leading term in the long-range expansion of the correlation function is calculated explicitly, 
and its gauge independence is proved. The leading contribution turns out to be of zero order in the 
Planck constant, and the relative value of the root mean square fluctuation of the Coulomb potential 
is found to be 1/V2, confirming the result obtained previously within the S-matrix approach. It 
is shown also that in the case of a macroscopic body, the h'^ part of the correlation function is 
- ; - suppressed by a factor 1 /N , where is the number of particles in the body. Relation of the 

^'bj^ I obtained results to the problem of measurability of the electromagnetic field is mentioned. 
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[Jh ■ I. INTRODUCTION 
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I Fluctuations in the values of physical quantities is an inahenable trait of aU quantum phenomena, which finds its 
reflection in the probabihstic nature of the quantum-mechanical description of elementary particle kinematics as well 
as of the field-theoretic treatment of the fundamental interactions. Investigation of fluctuations constitutes an essential 
part in establishing quasi-classical conditions, and therefore in dealing with the issue of correspondence between given 
' quantum theory and its classical original. 
7-H Precise identification of the quasi-classical conditions is especially important in investigation of the measurement 

process. In this process, an essential role is played by the measuring device "recording" the results of an observation, 
[ and the question of primary importance is to what extent this device can be considered classically. Furthermore, 
I investigation of fiuctuations in the observable quantities themselves is a significant ingredient in the theoretical treat- 
ment of the measurement process, in particular, in verifying conformity of the formalism of quantum theory with the 
i principal realizability of measurements. 
^ ' In the latter aspect, the issue of quantum fluctuations of the electromagnetic field was considered in detail by 
Oh, Bohr and Rosenfeld 1]. In this classic paper, a comprehensive analysis of the measurability of electromagnetic field 
^ ' was given, and compatibility of the restrictions following from the formal relations for the field operators with the 
actual limitations inherent to any process of measurement was verified, refuting the earlier arguments by Landau and 
Peierls 0] against this compatibility. In a subsequent paper (see also 0), the authors argued that the problem 
of charge-current measurements can be reduced to that already solved for the electromagnetic field. Later, a similar 
^ investigation was carried out in quantum theory of gravitation by DeWitt |^ . 

' An important point underlying investigation of Refs. P, 13 is that the question of measurability of the electro- 
magnetic field, and hence of charges and currents, can be considered proceeding from the commutation relations for 
the operators describing free electromagnetic field. In particular, these relations were used to estimate characteristic 
value of the electromagnetic field fluctuations, denoted in in by 6. Defined as the root mean square fluctuation, S is 
evidently of the order 0{\^)} It was argued in Ref. that this estimate holds true even in the presence of sources, 
provided that the charge and current distributions representing these sources allow classical description (which was 
also confirmed later by detailed calculations in Refs. [ESIESl)- 
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As far as one is interested in the electromagnetic field produced by the test bodies, this result is certainly sufficient 
to justify the use of the relation 6 = 0{^/h) : Since the choice of the experimental setup is at disposal of the observer, 
the measuring device may always be assumed "sufficiently classical." However, the state of affairs is different when 
the given electromagnetic field to be measured is considered. Classical assumption about the field producing sources 
is irrelevant in this case, and the above estimate does not apply. The problem of fundamental importance, therefore, 
is to determine fluctuations of flelds produced by nonclassical sources. 



As was shown in Ref. Q , the value of & varies depending on the ratio of the space and time intervals characterizing the measurement 
process, but in any case it is proportional to Below, our main concern will be the dependence of & on the Planck constant, so we 
suppress all other factors on which 6 may depend. 
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On various occasions, this issue has been the subject of a number of investigations. Electromagnetic field fluctuations 
springing up as a response of elementary particles' vacua to external fields or nontrivial boundary conditions were 
studied recently in pil. Illl llil Il4l | where references to early works can be found. Fluctuations of t he g ravitational 
field, induced by vacuum fluctuations of the matter stress tensor were considered in Refs. [l4L [la lla IitI . 

It should be noted that despite extensive literature in the area, only vacua contributions of quantized matter fields 
to the fluctuations of electromagnetic and gravitational fields have been studied in detail. At the same time, it is effects 
produced by real matter that are most interesting from the point of view of the structure of elementary contributions 
to the field fluctuation. 

The purpose of this paper is to investigate fluctuations of the electromagnetic field produced by a single massive 
charged particle, taking full account of quantum properties of the source-field interaction. It will be shown that the 
structure of fluctuations in this case is quite different from that of vacuum-induced fluctuations, or that found in the 
case of classical source. In particular, the root mean square fluctuation of the field turns out to be of zero order in 
the Planck constant, 6 — 0{hP), rather than 0(\/ft). It seems that this remarkable fact has not been noticed earlier. 
For a given electric charge of the source consisting of many particles, the contribution turns out to be inversely 
proportional to the number of constituent particles, so the classical 0{^/h) estimate is recovered in the macroscopic 
limit. 

The paper is organized as follows. In Sec.^]a preliminary consideration of the problem is given, and the Schwinger- 
Keldysh closed time path formalism used throughout the work is briefly reviewed. General properties of the correlation 
function of electromagnetic field fluctuations arc considered in Sec. lIIII In investigation of quantum aspects of particle- 
field interaction it is convenient to isolate purely quantum-mechanical effects related to the particle's kinematics 
assuming it sufficiently heavy. As discussed in Sec. Illl Al this leads naturally to the long-range expansion of the 
correlation function, making it appropriate to use the terminology and general ideas of effective field theories |l9l I20I 
|^|22. Next, it is proved in Sec. IIII Bl that the leading term in the long-range expansion of the correlation function, 
which describes fluctuations of the static potential, is independent of the choice of gauge condition used to fix the 
gradient invariance. Finally, the correlation function is evaluated explicitly in Sec. IIVI Sec. Ivl contains discussion of 
the results obtained. 



II. PRELIMINARIES 



Let us consider a single particle with mass m and electric charge e. In classical theory, the electromagnetic field 
produced by such particle at rest is described by the Coulomb potential 

^o = £:, A = o. (1) 

Our aim is to determine properties of this potential in quantum domain. For this purpose, it is convenient to separate 
this field-theoretic problem from purely quantum mechanical issues related to the quantum kinematics of the particle. 
Namely, the particle's mass will be assumed sufficiently large to neglect the indeterminacy in the values of particle 
velocity and position, following from the Heisenberg principle, so as to still be able to use the term "particle at rest." 

In quantum theory, Eq. is reproduced by calculating the corresponding mean fields, (^0), (^)- In the functional 
integral formalism, these are given by 

(4) = / exp{z5} , (2) 



where $ collectively denotes the fundamental fields of the theory, $ = {Afj_,(j),(j>*}, the components 4>,(l)* describing 
the charged particle which for simplicity will be assumed scalar, D$ is the invariant integral measure, and S — S[^] 
the action functional of the system. Integration is carried over all field configurations satisfying 



A^^O, cjj^-^q^^, {rr~.{<j>*r, for i^Too, 

where the superscripts + and — denote the positive- and negative-frequency parts of the fields, respectively, and 0o 
describes the particle state. Assuming that the gradient invariance of the theory is fixed by the Lorentz condition 



G EE d^A^ = , 



(3) 
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the action takes the form 



1 

2^ 



Fan — daAy — di/Afj 



(4) 



where ^ is the Feynman weighting parameter. Introducing auxihary classical sources 3 = {J'^ ^ j* , j} for the fields 
$ = {Afi,, (j>, (/>*}, respectively, the right hand side of Eq. (O may be rewritten in the standard way 



(4) = 
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exp 
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+j*{x)D{x,y)j{y)]+i I d'^x U(l^o+fM 



d''y[U^ix)D^,ix,y)J''iy) 



(5) 



3=0 



where S^^'' denotes the free field part of the action, S""*[<i>] = S[^] — S'^2^[$], and D are the propagators of the 
electromagnetic and scalar field, respectively, defined by 



d'^2 



SAf'{x)6A''{z) 

^2^(2) 



S(l)*{x)S(j){z 



D^^{z,y) = -S^^S^'^Xx-y), 
D{z,y) = -S'^^\x~y). 



(6) 
(7) 



It is not difficult to verify that the result of evaluation of the right hand side of Eq. (jsj in the tree approximation is 
given exactly by Eq. The field fluctuation, however, cannot be determined just as directly. The point is that the 
formal expressions like {Aq{x)) are not well defined because of the singular behavior of the product Ao{x)Ao{y) as 
X ^ y. This is a well known problem, encountered already in the theory of free fields. As was emphasized in Ref. Ill, 
in any field measurement in a given spacetime point, one deals actually with the field averaged over a small but finite 
spacetime domain surrounding this point, so the physically sensible expression for the field operator is the following 



^ I dt I d^'x A^{x,t). 
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Respectively, the product of two fields in a given point is understood as the limit of 

1 



V 



dt I dt' / d^x / d^x' A^,{x,t)A^{x' ,t') 



(8) 



(9) 
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when the size of the domain tends to zero. Finally, the correlation function of the electromagnetic 4-potcntial in this 
domain is 



(10) 



The necessity of spacetime averaging of field operators, Eq. dictated by the physical measurement conditions, 
entails the following important complication in calculating their correlation functions. The standard Feynman rules 
for constructing matrix elements of a product of field operators, such as that in the right hand side of Eq. (O, give 
the in-out matrix element of the time ordered product of the operators, i.e., the Green function, rather than its 
in-in expectation value. This circumstance is not essential as long as one is interested in evaluating (under stationary 
external conditions) the mean fields. However, it does make a difference in calculating the correlation function whether 
or not the field operators are smeared over a finite spacetime region. As is well known, in order to find the expectation 
value of a product of operators taken in different spacetime points, the usual Feynman rules for constructing the matrix 
elements must be modified. According to the so-called closed time path formalism 23, 24] (for modern reviews of the 
method, see Refs. [l^H^), the in-in matrix element of the product A^^{x)A'' {x') can be written as 



C'"'{x,x') = (in|>(a;)i''(a;')|in) 



D$+ Atix)A''+{x') exp{iS[^+] - iS[^-]} , 



(11) 
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where the subscript + (— ) shows that the time argument of the integration variable runs from — oo to +00 (from +00 
to —00). Integration is over all fields satisfying 



^+^0, 0i->(/<J, i^nl^i^oF for t 



-00 . 



$+ = for t +00 . (12) 

It is seen that (in|^'^ (a::)^'^(a::')|in) is given by the ordinary functional integral but with the number of fields doubled, 
and unusual boundary conditions specified above. Accordingly, diagrammatics generated upon expanding this integral 
in powers of the coupling constant e consists of the following elements. There are four types of pairings for each field 
or (j), corresponding to the four different ways of placing two field operators on the two branches of the time path. 
They are conveniently combined into 2x2 matrices^ 

■/?r+(x,y) Z?f_(x,y)\ A(r>(a;)i''(y))o z(iny)>W>o " 



^ f D++{x,y) D+^{x,y) \ ^ f i{T${x)$Hy))o t{${y)^Hx))o 
^ '^^ yD^+ix,y) D^^ix,y) J \ ^{${x)$Hy))o i{mx)4>Hy))o 

where the operation of time ordering T (T) arranges the factors so that the time arguments decrease (increase) from 
left to right, and (•)o denotes vacuum averaging. The "propagators" S^jy, D satisfy the following matrix equations 

d^z&,4^,z)T>''^z,y) ^ -tS'^^S^^\x-y), &^^{x,y) = i jj^^^^^^ , (13) 

dS0(x,z)S)(z,y) = -e^W(x-y), ©(x, y) = i ^^^^ , (14) 

where e, i are 2x2 matrices with respect to indices +, — : 

1 \ ■ _ 1 
1 j ' -1 

As in the ordinary Feynman diagrammatics of the S-matrix theory, the propagators are contracted with the vertex 
factors generated by the interaction part of the action, S''"*[$], with subsequent summation over (+, — ) in the vertices, 
each "— " vertex coming with an extra factor (—1). This can be represented as the matrix multiplication of Tifj_,y,D 
with suitable matrix vertices. For instance, the Adcfxp* part of the action generates the matrix vertex 2J which in 
components has the form 



'''' ' '■'^ 5A^{x)5cj>{y)5^*{z) 



where the indices i, j, k take the values +, — , and Sijk is defined by = — s = 1 and zero otherwise. External 

(f> {4>*) line is represented in this notation by a column (row) 



00 



\ 



satisfying 



d-^z ©(x,z)r(z) = , y^^x rt(x)©(x,z) = (0,0). (15) 



Figure n depicts the tree diagrams contributing to the right hand side of Eq. Hll|l . The disconnected part shown 
in Fig. ^a) cancels in the expression for the correlation function [see Eq. pOfi ]. which is thus represented by the 
diagrams (b)-(h). 



^ Below, Gothic letters are used to distinguish quantities representing columns, matrices etc. with respect to indices +, — 
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Explicit expressions for various pairings of the photon and scalar fields 



(27r)4 to2 - fc2 - iO ' — ^ '^^ J (27r)4 /fc2 „ ^2 _ jo ' 
D^+ix,y) = zj ^0ik'^)S{e-m^)e-^'^(--y\ D+^{x,y) ^ D^+{y,x) . (16) 

The photon propagator is written here in the Feynman gauge ^ = 1 which is most convenient in actual calculations. 
The question of how the choice of gauge condition affects the correlation function is considered in Sec. IIII Bl 

III. LONG-RANGE PROPERTIES OF CORRELATION FUNCTION 

Before we proceed to calculation of the correlation function, we shall examine its general properties in more detail. 
Namely, the structure of the long-range expansion of the correlation function, and the question of its gauge dependence 
will be considered. 

A. Correlation function in the long-range limit 

The mean electromagnetic field produced by a massive charged particle is a function of five dimensional parameters - 
the fundamental constants ft, c, the charge e and mass m of the particle, and the distance between the particle and the 
point of observation, r. Of these only two independent dimensionless combinations can be constructed - the constant 

/he playing the role of the expansion parameter of perturbation theory, and the ratio Ic/r = >c , where Ic = h/mc is 
the Compton length of the particle. As we have mentioned above, the particle is assumed sufficiently heavy so as to 
neglect effects related to the particle kinematics. This means that the quantum fluctuations are investigated in the 
limit X ^ Q. For fixed particle's mass this implies large values of r. In other words, the relevant information about 
field correlations is contained in the long-range behavior of the quantity C^^. 

To extract this information we note, first of all, that in the long-range limit, the value of C^j/ is independent of 
the choice of spacetime domain used in the definition of physical electromagnetic field operators, Eq. ((Sj). Indeed, 
in any case the size of this domain must be small in comparison with the characteristic length at which the mean 
field changes significantly. In the case considered, this means that V ^ . To the leading order of the long-range 
expansion, therefore, the quantity {A^{x)A^[x')) appearing in the right hand side of Eq. UlUI) can be considered 
constant within the domain. However, one cannot set a; = a;' in this expression directly. It is not difficult to see that 
the formal expression {A^{x)A^{x)) does not exist. Consider, for instance, the diagram ^b). It is proportional to 
the integral 



{x-x',p)= [ d^fcg(fc°)^(fc^)e"-(--') ^^^f + ^"^^"^s^r ^ y 

J [[k + q)'' — m^\[k — p)'^ 



where is the 4-momentum of the scalar particle, and the momentum transfer. For small but nonzero {x — x') this 
integral is effectively cut-off at large fc's by the oscillating exponent, but for x — x' it is divergent. This divergence 
arises from integration over large values of virtual photon momenta, and therefore has nothing to do with the long- 
range behavior of the correlation function, because this behavior is determined by the low-energy properties of the 
theory. Evidently, the singularity of I^,y{x — x',p) for x ^ x' is not worse than ln(a;^ — x'^^)"^ /{x^^ — x'^^Y. Therefore, 
{m\A^{x)A^{x')\\Ti) is integrable, and 6^,^ given by Eqs. (|10|l is well defined. 

Our aim below will be to show that this singularity can be consistently isolated and removed from the expression 
for If_t^{x — x',p) and similar integrals for the rest of diagrams in Fig.^ without changing the long-range properties 
of the correlation function. After this removal, it is safe to set a; = cc' in the finite remainder, and to consider C^^ as 
a function of the single variable - the distance r. An essential point of this procedure is that the singularity turns out 
to be local, and hence does not interfere with terms describing the long-range behavior, which guaranties unambiguity 
of the whole procedure. 

It should be emphasized that in contrast to what takes place in the scattering theory, the ultraviolet divergences 
appearing in the course of calculation of the in-in matrix elements in the coincidence limit are generally non-polynomial 
with respect to the momentum transfer. The reason for this is the different analytic structure of various elements 
in the matrix propagators S^iiy, S), which spoils the simple ultraviolet properties exhibited by the ordinary Feynman 
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amplitudes.'^ Take the above integral as an example. Because of the delta function in the integrand, differentiation of 
Ifii^ix — x' ,p) with respect to the momentum transfer does not remove the ultraviolet divergence of /^i^(0,p). What 
makes it all the more interesting is the result obtained in Sec. IIVI below, that the non-polynomial parts of divergent 
contributions eventually cancel each other, and the overall divergence turns out to be completely local. 

Next, let us establish general form of the leading term in the long-range expansion of the correlation function. 
In momentum representation, an expression of lowest order in the momentum transfer with suitable dimension and 
Lorentz transformation properties is the following 

e("jb)^4^- (17) 



Thus, only 00-component of the correlation function survives in the long-range limit: 



Not all diagrams of Fig.^contain contributions of this type. It is not difhcult to identify those which do not. Consider, 
for instance, the diagram (h). It is proportional to the integral 



{k-pY ' 

which does not involve the particle mass at all. Taking into account that each external scalar line gives rise to the 
factor (2£q)^^/2^ where Sq = \/ m? + q"^ « m , we see that the contribution of the diagram (h) is proportional to 1/m. 
The same is true of all other diagrams without internal scalar lines. As to diagrams involving such lines, it will be 
shown in Sec. llVl bv direct calculation that they do contain contributions of the type Eq. H17|l . But prior to this the 
question of their dependence on the gauge will be considered. 



B. Gauge independence of the leading contribution 

The value of 0-component of the momentum transfer is fixed by the mass shell condition for the scalar particle: 

pO = + {p + qY - v/m2 + = O f 



It follows from this estimate that is to be set zero in the long-range limit. This implies that Cj2 represents 
fluctuations in a quantity of direct physical meaning ~ the static potential energy of interacting particles. As such it 
is expected to be independent of the gauge condition used to fix the gradient invariance. More precisely, the issue of 
gauge dependence in the present case consists in the following. Although the gauge condition © fixes the form of 
the classical solution (QJ unambiguously, the inverse is not true: There are infinitely many ways of constructing the 
gauge fixed action, which lead to the same expression for the static potential. We have to verify that this freedom in 
the choice of the gauge condition does not affect the value of C^'J] . 

The gauge independence of the leading contribution will be demonstrated below in the most important particular 
case when changes of the gauge-fixed action are induced by variations of the Feynman weighting parameter ^. It is 
not difficult to see that these variations do not alter the value of classical potential. Indeed, the latter satisfies 

(IS, 

where Sera is the gauge-fixed action for electromagnetic field. 



^ As is well known, the proof of locality of the S-matrix divergences relies substantially on the causality of the pole structure of Feynman 
propagators. This property allows Wick rotation of the energy contours, thus revealing the essentially Euclidean nature of the ultraviolet 
divergences. 
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and J'^ the electromagnetic current, 

Acting on Eq. (fTS|l with 9^, using the current conservation 9^ = and the identity d^duF^^" — 0, we obtain 

which iniphes that G = on the classical solution, and that ^ effectively falls off from the classical equation (I18f) . Let 

us now prove that C^°j remains unchanged under variations of ^ too. We shall proceed along the lines of Ref. 
where ^-independence of hP part of the one-loop correction to Newtonian potential was proved within the S-matrix 
approach. 

As we saw in Sec. IIII Al the only diagrams contributing in the long-range limit are those containing internal scalar 
lines. We will show presently that the ^-dependent part of these diagrams can be reduced to the form without such 
lines. First, it follows form Eq. that the ^-derivative of the electromagnetic propagator satisfies 



On the other hand, acting on Eq. (|13|) by 9^ one finds 

jIU^d:^,^{x, y) = -EdlS^^^ {x~y). (20) 

This equation can be resolved with respect to d'^'Ziua- The reciprocal of the operator /□, appropriate to the boundary 
conditions \i2l is (— S'^), hence 

d'i^,o.{x,y)=^dl^\x,y). (21) 

Taking into account that 

substituting this into Eq. 119|l . and using Eq. (|21|) . we get 

^^M^l^ = -i| d^z S),„(x,z)/a,"9^S)°(.,y), 

or 

^^^g^""'^^ = d^dlaix, y) , a(x, y) ^ J d'z ^^(x, z)/DO(z, y) . (22) 

Thus, assuming the matrix indices and spacetime coordinates referring to the point of observation fixed, we see that 
the ^-dependent part of the matrix photon propagator is contracted with the matrix vertex 03^ through a gradient 
term 

dM^)-{f^b'it^)) ' K:^)=5>(x,y). 

Integrating by parts the spacetime derivative may be rendered to act on the vertex. 

Next, we use gauge invariance of the action S'o[$] to transform bd^JB^. This invariance is expressed by the identity 

Differentiating Eq. (|23|l twice with respect to </>, 4>* , and setting $ gives 



d: 



^'5A^{x)5<jy{y)5<P*{z) 



- . ,ie5^^\x - y) ~ , ° , A e5^^\x-z). (24) 
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The matrix vertex is obtained by multiplying S^So/5(j)S4>*SAf^ by the matrix s^fe. It follows from Eq. 124() that 
the combination b9^2J^ may be written as 



:2 o(2) ^ 



b9^2J'^(x, y, z) = ie [s{z, x)S^*^ {x - y) ~ d{x, y)^^ (x - z) } , 



(25) 



where 



p c.(2) 



Sr{x)6cf>{y) ' 



6+ 

' ^ -6- 



Finally, contracting bdf^'V^{x,y, z) with the matrix 'D{y,y') and the vector r(z) = (0o(2:), ^ol^)) , and using 
Eqs. (O, (Unjl gives 



ief(a;) d'^y^{x,y)D{y,y') = ie^oM^''^ (x - y') • 



We see that upon extracting the ^-dependent part of the diagrams (b)-(f) in Fig.^ the scalar propagator is cancelled 
by the vertex factor. Using the reasoning of Sec. IIII Al the resulting diagrams give rise to terms ~ 1/m. Thus, 
^-independence of the correlation function in the long-range limit is proved. 

IV. EVALUATION OF THE LEADING CONTRIBUTION 

Let us proceed to calculation of the leading contribution to the correlation function. It is contained in the sum of 
diagrams (b)-(f) in Fig.^ which has the symbolic form 



where the superscript "tr" means transposition of the indices and spacetime arguments referring to the points of 
observation: + ^—,x<^-x' [the transposed contribution is represented by the diagrams collected in part (f) 

of Fig.^. As it follows from the considerations of Sec. IIII c|S can be expressed through C^i, as 



Written longhand, If^^, reads 



el°j = hm I T-^-T^ 



d xd x' C^^{x, x') 



{V.T) 



(26) 



/^^(x,x') = ?e2 / / d^zd^z'l + 





{x 


^) 


r(^)i 


D+ 




z') dl q 


Hz') 






x') 




[x 


^) 


0*(;^)i 


D+ 


-iz, 


z') dl c, 


h{z') 


D'L_ 




,x') 




[x 


^) 


r(^) ^ 


D_ 


+ {z 


z') 5 <; 


b{z') 


Dl_ 




,x') 


D\_ 


[x 


^) 






-iz 


z') 5 ^ 


6(z') 






,x') 



(27) 
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where 

Contribution of the third term in the right hand side of Eq. (|27|l is zero identicahy. Indeed, using Eq. and 
performing spacetime integrations we see that the three hnes coming, say, into z-vertex are ah on the mass shell, 
which is inconsistent with the momentum conservation in the vertex. The remaining terms in Eq. (|27() take the form 

r , f^ 2 ff d^P a*{q)a{q + p) ~ . . 

= e yy ^^_=_e po=e,+p-£,, (28) 

where 

V(P.9)=-* / e-'=(--')(2g^ + fc^)(2g, + fc,+p,) 

x{ D\^{k)D^^{q^k)D\_(k-v) 
-Dl+ik)D+^{q + k)D'L_{k ~ p) 

+D'l_{k)D^^{q + k)D"__{k-p) }, (29) 
and a{q) is the Fourier transform of the particle wave function, normalized by 

'a(q)r = 1- 



(2^) 



The function a{q) is generally of the form 

where Xq is the mean particle position, and b{q) is such that 



J dPq b{q)e"i'' = (30) 



for X outside of some finite region W around x = 0. 

Upon extraction of the leading contribution these expressions considerably simplify. First of all, one has e^+p ~ 
£q « m, and hence, po ~ 0. Next, b{q + p) may be substituted by b{q) : This implies that we disregard spatial 
spreading of the wave packet, neglecting the multipole moments of the charge distribution. Taking into account the 
normalization condition, Eq. H28|) thus becomes 

= (03(^I^(9)^^"^^"^"^V(^^,9), Po = 0. (31) 

The leading term in /^^ (p, q) has the form [Cf . Eq. (fTTjl ] 

ii°)(p,g)^44=- (32) 

This singular at p ^ contribution comes from integration over small k in Eq. (|29|l . Therefore, to the leading order, 
the momenta k,p in the vertex factors can be neglected in comparison with q. Thus, 

V(p,g) = ^Aq^q.i J ^e-'^(---'){ Dl4k)D++{q + k)Dl_{k-p) 

- Dl+{k)D+^{q + k)D^__{k-p) 

+ Dl_{k)D^^{q + k)D°__{k-p)Y (33) 
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Furthermore, it is convenient to combine various terms in this expression with the corresponding terms in the trans- 
posed contribution. Noting that the right hand side of Eq. \6'6\i is exphcitly symmetric in /i, z/, and that the variables 
x,x' in the exponent can be freely interchanged because they appear symmetrically in Eq. f^ . we may write 



X { [Dl4k)D++{q + k)Dl_{k -p) + D°__{k)D--{q + k)D''_+{k - p)] 
+ [Dl_{k)D^^{q + k)D"__{k -p) + D\{k)D++{q + k)D\4k^p)\ 
- [D\4k)D+^{q + k)D'l_{k-p) + D'l_{k)D^+(q + k)D\+{k-p)\]. (34) 

With the help of the relation 

£'__+£»++ = £'+_+£»_+, (35) 

which is a consequence of the identity 

T4>{x)^{y) + f4>{x)4>{y) - 0(x)0(y) - 0(y)(^(a;) = , 

the first term in the integrand can be transformed as 

D%4k)D++{q + k)D%_{k-p) + D"__{k)D--{q + k)D°_+{k - p) 
= Dl^{k)D++{q + k)D°{k -p) + D°{k)D__{q + k)D''_^{k ~ p) , 

where 

D{k) EE D+^{k) + D^+{k) = 2m5{k^ - m^) , D°{k) = D{k)\m=o ■ 

Here we used the already mentioned fact that D^_{k)D |-(fc + g) = for q on the mass shell. Analogously, the 

second term becomes 

D^{k)D^^{q + k)D"__{k-p) - D\(k)D^^{q + k)D°{k-p) . 

Changing the integration variables k—^k + p^q^q—p, and then p —p, noting that the leading term is even 
in the momentum transfer [see Eq. H32|l ]. and that p(x — x') in the exponent can be omitted in the coincidence limit 
[x x'), the sum of the two terms takes the form 

D\k) [D++{q + k)Dl^{k -p) + D^^iq + k)D^__{k^p)] . 

Similar transformations of the rest of the integrand yield 

D%ik)D+^{q + k)D'l_{k^p) + D'l_{k)D^+{q + k)D\^{k - p) 

^ -i [D\+{k)D{q + k)Dl+{k -p) + D'L_{k)D{q + k)D'L_{k - p)] . (36) 

Substituting these expressions into Eq. H34(l and using Eq. H16|l gives 

d'^k 
(2^ 

^ R,, , S[{q + kf-m\ , 

' [(/fc-p)2+i0][(g + fc)2-m2 + z0] ^ [k^ + iQ\[{k-pY ' ^ ' 



ct^k 



As was discussed in Sec. lIII Al the exponent in the integrand in Eq. (|37|l plays the role of an ultraviolet cutoff, ensuring 
convergence of the integral at large fc's. On the other hand, the leading contribution l|32|) is determined by integrating 
over fc ~ p where it is safe to take the limit x — > x' . Since (x — a;') is eventually set equal to zero, one can further 
simplify the k integral by using the dimensional regulator instead of the oscillating exponent. Namely, introducing 
the dimensional regularization of the k integral, one may set x = x' afterwards to obtain 



"(27rp" I I 

^ (38) 



(27r)3 \ [{k - pY + iO] [{q + ky -~m? + iO] 
<5[(g + fc)2 - m2] 



[k"^ + iQ\[{k - pY + iQ] 
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where fj, is an arbitrary mass parameter, and e = 4 — d, d being the dimensionahty of spacetime. 
Next, going over to the a-representation, the first term in the integrand may be parameterized as 



Substituting this into Eq. I|38(l . and using the formulas 



[{k - J3)2 + iO] [{q + kf -m^ + iO] (p2 - 2kp + iO) {2kq + iO) 



1 / TT 



i \\a\ 



2-e/2 







exp 


(- 






7re\ 


H 


27/ 


a 



one finds 



K{p) = fi' 



d^-^k 



(27r)3 [{k - p)2 + iO] [{q + fc)2 - 



(2^) 



/i%-^/2 



■ COS 



dxdydz ^— ^ 

(?) 



exp 



r 1- 



dydz e 



lyp- 



-{yp- zqf 



2 \ {yp- zqfl"'"^ 



exp 



-{yp~ zqf 



87r2 ™V4y^V^ 27 7^0 " {yp~zq) 

Changing the integration variable z ^ yz, and taking into account that q'^ — m? , qp — —p'^/2 yields 

|z2_(l + z)ap/2 



(p) = — — ^ cos 



-e/2 



(t) 



r(l-- / dye'^^'y-i / dz^ ^ 



^2 



cos 



|z2 - (1 + z)a 

87r2m2 \-p'^ J V 4 y " V" 2)^"'"']^ "'~ z2 - (1 + z)a 
where a = —p^/m? . Similar manipulations with the second term in Eq. H^^8(l give 



e/2 



d^-'k (5[(g + fc)2-m2 



dxdydz 



(27r)3 [fc2+iO][(fc-p)2+iO] (27r)4 777o 

d^'^'fe |^e"(*:^+2fc</) _|_ g-ix(fc^+2fcg)~j giy(p^-2fc(g+p) + i0)gi2:(-2fc9+i0) 
J^" j j j ' ' ^""^ 



(2^) 



exp 



/TT 

dxdydz ( — 



„i7re/4 



exp 

-e/2 



167r2m2 (m) ( ^ 2 

-i7re/4 



X /y^ '^y'^^ \[{y + z + iy+yzaY+^/^ ' [(y + z- 1)2 + yza] 1+^/2 



Changing the integration variables y — ut, z = {1 — u)t yields 



Hp) 



(9^(-|) 



(39) 



(40) 



L "^Vo '^'^ U [(t+ 1)2 +W(l-U)t2a] 1+^/2 ' [(t-l)2+w(l-u)i2c,]l+e/2 

On the other hand, L{p) must be real because the poles of the functions £>++(&), D++(fc — p) actually do not 
contribute. Hence, 



Lip) 



e/2 



(£)'»(?)r(-|) 



dt 



t du 



7o [(t+l)2+ii(l-w)t2a]i+'^/2 
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Let us turn to investigation of singularities of the expressions obtained when e — > 0. Evidently, both K{p) and L{p) 
contain single poles 



Vi- 


f 4/a4 


- 1 


vi- 




- 1 


an 


-4/a + 


1 


v/n 


-4/a- 


1 



"''^P' = 8^^e Jo [l + uil~u)a] " An^m^e ^/TTIJ^ ""^T+t/T^ ' ^^^^ 

Note that L^w — — 2ifdiv Upon substitution into Eq. (|38|) the pole terms cancel each other. Thus, Cf^^{p,q) turns 
out to be finite in the limit e ^ 0. Taking into account also that divergences of the remaining two diagrams (g), (h) 
in Fig. Hare independent of the momentum transfer,^ we conclude that the singular part of the correlation function 
is completely local. 

It is worth of mentioning that although the quantity ii^div(p) is non-polynomial with respect to p^, signifying non- 
locality of the corresponding contribution to C^^(x, x'), it is analytic at p = 0, which implies that its Fourier transform 
is local to any finite order of the long-range expansion. Indeed, expansion of K^iv{p) around a = reads 

KMp) = -T7r^(^~^ + ---+^-^f^^'' + ---) ■ (42) 
lD7r^r7?/e \ 6 (2n-|-l)! J 

Contribution of such term to C^v{x, x') is proportional to 

The delta function arose here because we neglected spacial spreading of the wave packet. Otherwise, we would have 
obtained 

(tr '^^'^^'^'^ + p)e'^(-'-") (l ~ ^ + ■■■)- 1- x'-Xo^W, 

as a consequence of the condition H3U|) . In particular, applying this to the correlation fmiction, we see that its divergent 
part does not contribute outside of W. Thus, the two-point correlation function has a well defined coincidence limit 
everywhere except the region of particle localization. 

Turning to calculation of the finite part of K{p), we subtract the divergence from the right hand side of Eq. H4U|I. 
and set e = afterwards: 



iffi„(p) = lim[if(p)-ifdivb)] 



1 



j'tt + In TT + 7 + 2 In 



2 



f 

Jo 



dz 



2^ - (1 z)a 



f 

Jo 



ln|z2- (l + z)a| I , , 

— (1 + z)a 



where 7 is the Euler constant. The leading contribution is contained in the last integral. Extracting it with the help 
of Eq. (|A3|) of the Appendix, we find 

K^niP) = -r^^- (44) 



The corresponding integrals do not involve dimensional parameters other than p, and therefore are functions of only. Hence, on 
dimensional grounds, both diagrams are proportional to 

(t-»)(^)"' = t-' + tK^)*°'* 

where ci^2 are some finite constants. 
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As to L{p), it does not contain the root singularity. Indeed, 

7r-'/2 f fiy (T^<^\ r(l + e/2) 



^(0) = ^V^f-) cos(^) 
^ ' STT^m^ \mJ V 4 / 



and therefore Lfin{p) = lim[L(p) — Ldiv(p)] is finite at p = 0. It is not difiicult to verify that Lfin(p) is in fact analytic 

e— ►O 

" , ■ : . Thus, substituting Eq. 144|l into Eqs. (|38(l . H31(l . and then into Eq. I|26() . and using the formula 

d^p e^P^ 1 



(27r)3 IpI 27r2a;2 ' 

we finally arrive at the following expression for the leading long-range contribution to the correlation function 

32^r2 



eo?(^) = T;7^, r^\x-Xo\, (45) 



all other components of cfS vanishing. This result coincides with that obtained by the author in the framework of 
the S-matrix approach 28] . The root mean square fiuctuation of the Coulomb potential turns out to be 



(A2(r)> = -=^. (46) 
V32 7rr 



Note also that the relative value of the fiuctuation is 1/ ^/2. It is interesting to compare this value with that obtained 
for vacuum fluctuations. As was shown in Ref. the latter is equal to V2 (this is the square root of the relative 
variance used in Ref. 0|). 

We can now ask for conditions to be imposed on a system in order to allow classical consideration of its electro- 
magnetic interactions. Such a condition can easily be found out by examining dependence of the /i" contribution on 
the number of field producing particles. Let us consider a body with total electric charge Q and mass M, consisting 
of a large number N = Q/q of elementary particles with charge q and mass m, assumed identical for simplicity. Then 
it is readily seen that the diagrams (b)-(f) in Fig. ^ are proportional to N ■ q'^ — Q'^/N because they have only two 
external matter lines. We conclude that the contribution to the correlation function, given for an elementary par- 
ticle by Eq. H45|l . turns into zero in the macroscopic limit N —> oo. At the same time, the next term in the long-range 
expansion of the correlation function has the form 

Hence, for a fixed mass M of the multi-particle body, it is independent of N, 

(M/N) cr^ M cr^ ' 

thus recovering the classical estimate S = 0{Vh) in the macroscopic limit. 



V. DISCUSSION AND CONCLUSIONS 



The main result of the present work is the expression H46() describing quantum fiuctuations of the electrostatic 
potential of a charged particle. It represents the leading contribution in the long-range expansion, and is valid at 
distances much larger than the Compton length of the particle. Remarkably, the fluctuation turns out be zero order 
in the Planck constant. This result is obtained by evaluating two-point correlation function of the electromagnetic 4- 
potential in the coincidence limit. We have shown that despite non-locality of divergences arising in various diagrams 
in this limit, the total divergent part of the in-in matrix element is local, and hence does not contribute outside the 
region where the particle is localized. 

Perhaps, it is worth to stress once more that the issue of locality of divergences is only a technical aspect of our 
considerations. This locality does make the structure of the long-range expansion transparent and comparatively 
simple. However, even if the divergence were nonlocal this would not present a principal difficulty. A physically 
sensible definition of an observable quantity always includes averaging over a finite spacetime domain, while the 
singularity of the two-point function, occurring in the coincidence limit, is integrable (see Sec. IIII A|l . The only 
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problem with the nonlocal divergence would be impossibility to take the limit of vanishing size of the spacetime 
domain [Cf. Eq. fl^]. 

Turning back to the problem of measurability of electromagnetic field, which served as the starting point of our 
investigation, we may conclude that the proof of principal rcalizability of measurements, given in Refs. P, is 
incomplete, since it relies heavily on the 0{Vh) estimate for the electromagnetic field fluctuations, which is inapplicable 
to the fields of elementary particles. The latter case therefore requires further investigation. 

As we saw in the end of Sec. IIVI the contribution to the correlation function of the electromagnetic field of a 
macroscopic body is suppressed by the factor 1/A^. This is in agreement with the necessary condition for neglecting the 
ff' contribution, formulated in Ref. 29] in connection with the loop corrections to the post-Newtonian expansion in 
general relativity. In the present context, it coincides with the well-known criterion for neglecting quantum fluctuations 
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APPENDIX A 



In the course of separation of the finite part in the right hand side of Eq. 14U|I we encountered the integral 
^ r-^^ln^^_(l±iH ^ In|(.-P.)(.-P.)| 



z2-(l + z)a Jq {z-pi){z-p2) 

a ± + 4a 

P1.2 = ^ . 

To evaluate this integral, let us consider the integral 

Jc [W - pi)(W - P2) 

taken over the contour C in the complex w plane, shown in Fig. [3 A is zero identically. On the other hand, one has 
for sufficiently small positive a 

. fP'-\ \n[(z-p,)(z-p2)] . Hp2-Pi) iiT\na + 7Ty2 

A — dZ ; ^ ITT- 



>-oo {Z-Pl){z-P2) P1-P2 P1-P2 

rP2-a \n\(z - pi){jp2 - z)] - i-K . In(p2 - Pi) - iTT iTrlncT + 7r2/2 
-I- / dz 7 m 

Jp^+a [Z-Pl){z-P2) P2-P1 P2 - Pi 

^ r H{z-p,){z^p2)]-2tn 

Jp2+a {Z-Pl){z-P2) 

Changing the integration variable z — > (pi + P2) — z in the first integral, and rearranging yields in the limit cr ^ 

A . i d. + ^ . (A2) 

2 Jo \z-pi){z-p2) P2-P\ 

The root singularity is contained in the second term, because 

^^+^^ \n\{p2-z){z-p^)\ 

dz — ^7 > — ma tor a — > U. 

\z - V\)K.z ~ P2) 

Thus, 

(A3) 



4Va ■ 
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It is interesting to note also that the part containing the root singularity [the second term in Eq. (|A2|) ] exactly 
coincides with that found using the S-matrix (Cf. Eq. (B2) in Ref. j2q'). 
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FIG. 2: Contour of integration in Eq. lAH . 



